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Re-analysed JADE data were used to determine as at y/s = 14-44 GeV on the basis of resummed calculations 
for event shapes and hadronisation models tuned to LEP data. The combined result is Qs(M z o) = 0.1194+o'oogg 
which is consistent with the world average. Event shapes have also been used to test power corrections based on 
an analytical model and to verify the gauge structure of QCD. The only non-perturbative parameter cto of the 
model was measured to ao(2 GeV) = 0.503 t o 045 an d is found to be universal within the total errors. 



1. INTRODUCTION 

The re-analysis of e + e~ annihilation data col- 
lected with the JADE detector at the PETRA col- 
lider (1978-1986) has been shown to be a valuable 
effort Q^,&H] since the characteristic energy evo- 
lution of Quantum Chromodynamics (QCD) be- 
comes more manifest towards decreasing centre- 
of-mass energies y/s. Recently, data at energies 
down to s/s = 14 GeV could be employed in state- 
of-the-art QCD studies due to the successful res- 
urrection of the original JADE software. 

Since the PETRA shutdown, significant 
progress has been made in the theoretical calcu- 
lations of event shape observables. In the follow- 
ing, we present an as analysis at y/s = 14-44 GeV 
based on the most complete perturbative calcu- 
lations for event shapes m available so far. We 
included recently analysed data in the energy re- 
gion = 14-22 GeV for which the calculations 
are applied for the first time. Furthermore, power 
corrections based on an analytical model by Dok- 
shitzer, Marchesini, and Webber [Q (DMW) were 
investigated as a promising approach to describe 
non-perturbative effects in event shapes. Besides 
as, the model depends only on one additional free 
parameter. Also the consistency of power correc- 
tions with the gauge structure of QCD was tested. 

2. EVENT SHAPES 

From multihadronic data samples, the distri- 
butions of thrust (1 — T), heavy jet mass (Mh), 
total and wide jet broadening (£?t and -Bw), C pa- 



rameter and the differential 2-jet rate j/23 in the 
Durham scheme are calculated (cf. Q ) . The data 
are corrected for the limited acceptance and res- 
olution of the detector and for initial state pho- 
ton radiation. Since electroweak decays of the 
heavy b-hadrons fake hard gluon radiation par- 
ticularly at yfs < 22 GeV, we take the contribu- 
tion e + e _ — > bb as an additional background to 
be subtracted from the distributions. 

The data have been used to assess the perfor- 
mance of various QCD event generators tuned 
to LEP data at «/s = M z o (cf. g). The par- 
ton shower and string fragmentation model im- 
plemented in Pythia/Jetset is found to be 
well capable of describing event shapes down 
to 14 GeV. The quality of the models Ariadne 
(colour dipole scheme) and Herwig (cluster frag- 
mentation), however, is more moderate, and 
Cojets (independent fragmentation) is clearly 
disfavoured. Obviously, the model parameters of 
these generators need a re-tune at lower s/s. 

3. DETERMINATION OF a s 

The determination of ag is based on a com- 
bination of an exact QCD matrix element cal- 
culation 0{ol\) § intended to describe the 3- 
jet region of phase space and a next-to-leading- 
logarithmic approximation (NLLA) || valid in 
the 2-jet region where multiple radiation of soft 
and collinear gluons from a system of two hard 
back-to-back partons dominate. We perform \ 2 - 
fits of the theoretical predictions corrected for 
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Table 1 

Preliminary as results from JADE. 



a. 



[GeV 



as (y/s) fit+exp. had. hi. ord. 



tot. 



14.0 0.1704 

22.0 0.1513 

34.8 0.1431 

38.3 0.1397 

43.8 0.1306 



±0.0051 
±0.0043 
±0.0019 
±0.0040 
±0.0037 



+0.0141 
-0.0136 

±0.0101 
±0.0073 
±0.0054 
±0.0056 



+0.0143 
-0.0091 
+0.0101 
-0.0065 
+0.0091 
-0.0060 
+0.0084 
-0.0056 
+0.0068 
-0.0044 



+0.0206 
-0.0171 
+0.0144 
-0.0121 
+0.0118 
-0.0096 
+0.0108 
-0.0087 
+0.0096 
-0.0080 



hadronisation effects. For the main results, we 
use the ln(i?)-scheme [|| for the perturbative pre- 
diction with the renormalisation scale factor = 
n/y/s = 1 and Pythia for the estimation of non- 
perturbative contributions. We generally observe 
stable fits and good agreement with the data at 
each y/s with y 2 /d.o.f. ~ 0.2-2.0. In case of B\\r, 
a significant excess of the theory over the data in 
the 3-jet region of the distributions is present. 

Experimental errors are under control for all 
data samples. Expectedly, hadronisation uncer- 
tainties increase rapidly towards 14 GeV. The in- 
dividual results agree with each other within 1- 
2 standard deviations of the fit and experimen- 
tal errors. For each y/s, the as values from the 
six observables are combined using the weighted 
mean method of Ref.0 (Tab.|). The total errors 
are dominated by higher order uncertainties. At 
14 and 22 GeV, hadronisation uncertainties are of 
the same order as the QCD scale ambiguities. 

The as results obtained here and in similar 
analyses at higher energies based on resummcd 
event shapes (Fig. [l]) agree well with the QCD 
expectation for the running coupling ||. A y 2 - 
fit taking statistical and experimental errors into 
account yields a s (Af z o) = 0.1213 ± 0.0006 with 
X 2 /d.o.f. = 8.3/11. Even considering the total 
errors, the unphysical hypothesis as = const, is 
disfavoured by a fit probability of « 10 -5 . 

4. TEST OF POWER CORRECTIONS 

The DMW model §] describes non- 
perturbative effects to event shapes as contribu- 
tions from gluon radiation at low energy scales, 
assuming that the physical strong coupling as (/i) 
remains finite in the energy region around the 
Landau pole where simple perturbative evolution 



s 
0.2 



0.15 



0.1 



J JADE (preliminary) 

J other experiments 

(e e jet rates & event shapes) 




QCD: a s (M z ) = 0.1184 ± 0.0031 
S. Bethke, J.Phys.G26 (2000) R27 



10 



50 



100 



200 

VsfGeV] 



Figure 1. a>s(y/s) derived using 0(a§)+NLLA pre- 
dictions for event shapes in e + e~ annihilation. 



of as breaks down. This leads to the introduction 
of a parameter ao(^i) = l///i / W d/u, as (/i) that 
absorbs all non-perturbative details of ag(fx) up 
to an arbitrary infrared matching scale fi\. The 
principle structure of power corrections is a shift 
VDjr of the perturbative spectrum away from the 
2-jet region, with V oc ao/y/s and Djr depending 
on T . In case of Bt and B\v, the shift is super- 
imposed by a squeeze Dp oc In 1 j T . For mean 
values one obtains also an additive correction. 
The perturbative part is C(ag)±NLLA for the 
distributions and 0(a%) for the means. 

The prediction, in particular the universality 
of ao(/j.i), has been tested by global fits to the 
hadron level data from this analysis and from 
other experiments e.g. at LEP/SLC up to y/s = 
189 GeV (cf. @), with a s (M z o) and a (2 GeV) as 
only free parameters. The major features of the 
distributions are reproduced well, thus support- 
ing the 1/y/s evolution of power corrections in 
event shapes. However, we observe discrepancies 
for the distributions of the less inclusive quanti- 
ties Mh and Bw m particular at y/s < M z o . 

For y23 the leading power correction is known 
to be quadratic in 1/ y/s. This expectation has 
been verified by means of the new JADE distri- 
butions at y/s = 14 and 22 GeV using a simple 
additive power correction ansatz. 

As shown in Fig. ||, there is a reasonable agree- 
ment between the individual results within the 
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Figure 2. cts(M z o) and ao(2GeV) and one standard 
deviation errors of DMW-fits to distributions a) and 
means b) . The hatched band represents the combined 
as{M z a) derived from the "conventional" analysis. 

total uncertainties. However, the ato results from 
the Mh and % distributions are large compared 
to the results from the other observables. This ob- 
servation may be related to the non-inclusiveness 
of these variables. The as values from power 
corrections to the distributions are systematically 
smaller than the results based on MC corrections. 
This is due to the different amounts of squeeze 
of the perturbative spectrum predicted by both 
types of model, particularly in case of £?w 

Combining the results for the mean values and 
the distributions taking correlations between the 
systematic errors into account yields as(M%a) = 
0.1175±°;°°3i an d a (2GeV) = 0.503±°;°45. The 
scatter of the ao values is mostly covered by the 
theoretical uncertainty of the Milan factor Q. 

5. STUDY OF QCD COLOUR FACTORS 

The DMW ansatz has been exploited to ex- 
tract the QCD colour factors Ca, Cp, and nf jjj. 
Various global fits to the event shape spectra try- 
ing alternative sets of the free model parameters 
support the SU(3) symmetry group. The most 
stable and precise measurements are provided by 
1 — T and C. Combining the corresponding re- 
sults for these variables with ao(2GeV) and rif 
fixed, and as(M z o), Cp, and Ca free, one finds 
C F = 1.29±0.18 and C A = 2.84±0.24. This is in 
good agreement with the QCD expectation while 
some other gauge symmetry groups are excluded. 



6. CONCLUSIONS 

Resummed QCD theory combined with LEP 
tuned hadronisation models fits event shape data 
well down to y/s = 14GeV and allow consis- 
tent determinations of cvg- The combined re- 
sult evolved to the Z° mass scale is as(M z a) = 

0. 1194.0.0068 w hich is substantially more precise 
than former PETRA measurements and also in 
good agreement with the world average value || . 

Power corrections cx 1/ yfs generally reproduce 
the overall event shape spectra, except for the dis- 
tributions of the less inclusive variables (Mr and 
£>w) at ^/s < AT z o. The results for a support 
the DMW prediction of universality within 25%. 
Using power corrections, the gauge structure of 
QCD has been verified with uncertainties com- 
petitive e.g. with traditional 4-jet angular corre- 
lation analyses. Potential biases from hadronisa- 
tion models are reduced within this approach. 

Thus, exploiting JADE data significantly im- 
proves the verification of QCD on the basis of 
e + e _ annihilation. 
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